arXiv: 1509.03963v2 [quant-ph] 14 Dec 2015 


NMR Investigation of the Quantum Pigeonhole Effect 


Anjusha V. S., Swathi S. Hegde, and T. S. Mahestj^ 

Department of Physies and NMR Researeh Center, 

Indian Institute of Seienee Edueation and Researeh, Pune 4^^008, India 

NMR quantum simulators have been used for studying various quantum phenomena. Here, using 
a four-qubit NMR quantum simulator, we investigate the recently postulated quantum pigeonhole 
effect. In this phenomenon, a set of three particles in a two-path interferometer often appears to 
be in such a superposition that no two particles can be assigned a single path, thus exhibiting the 
non-classical behavior. In our experiments, quantum pigeons are emulated by three nuclear qubits 
whose states are probed jointly and noninvasively by an ancillary spin. The experimental results 
are in good agreement with quantum theoretical predictions. 
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1. INTRODUCTION 

Quantum theory is well known for generating scenarios 
having no classical analog. A recent addition to this list 
is quantum pigeonhole effect (QPHE) [H |2]. In math¬ 
ematics, the pigeonhole principle states that if n items 
are put in m < n containers, then at least one container 
must have more than one item [3]. While the method 
of contradiction is sufficient for proving the pigeonhole 
principle for a finite number of objects, it calls for more 
advanced methods in the case of an infinite set of objects 
and containers. The generalized concept has several in¬ 
teresting consequences in mathematics mi, computer 
science mi, graph theory uni, and combinatorics mi. 

Recently, Aharonov and coworkers have theoretically 
illustrated certain quantum mechanical scenarios appear¬ 
ing to contradict the pigeonhole principle [T]. This phe¬ 
nomenon, known as QPHE has already raised consider¬ 
able interest. Eor example Yu and Oh demonstrated the 
emergence of QPHE from quantum contextuality [12]. 
Rae and Eorgan suggested that QPHE arises as a re¬ 
sult of interference between the wavefunctions of weakly- 
interacting particles m- In this work, we simulate 
QPHE using a four-qubit NMR quantum simulator. 

Molecular nuclear spins in liquid ensembles controlled 
by NMR techniques have successfully been employed for 
simulating various quantum phenomena. They include 
simulating quantum harmonic oscillators El, quantum 
phase transitions m, ground state of hydrogen atom 
HE], probabilities of particle in a potential El, and quan¬ 
tum delayed choice experiment m- In this work we 
utilize a 4-qubit quantum register based on an ensem¬ 
ble of four mutually interacting spin-1/2 nuclei partially 
oriented in a liquid crystal matrix. While the rapid 
isotropic reorientation of molecules in a liquid erases 
all the anisotropic interactions, only weak indirect spin- 
spin interactions (J-couplings) survive resulting in long 
nonlocal quantum gates that are comparable to deco¬ 
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herence time-scales. On the other hand, in solid-state 
NMR systems, the strong intra- and intermolecular in¬ 
teractions inhibit high-fidelity quantum control of desired 
qubits. In this sense, the partially oriented spin systems 
with no translational order have a specific advantage, 
i.e., negligible intermolecular interactions and residual in¬ 
tramolecular anisotropic interactions just strong enough 
for quantum control m- Accordingly, such systems are 
increasingly being used for realizing quantum testbeds 

umiiiso]. 

In the following section we provide a brief theoretical 
description of QPHE. In section III, we shall describe 
its NMR quantum simulation and finally we conclude in 
section IV. 


2. THEORY 

Let us first consider a single quantum particle entering 
the Mach-Zehnder interferometer shown in Eig. [^(a). It 
involves two beam-splitters (BSI and BS2), a 90 degree 
phase-shifter (Z), and two detectors (DO and DI). BSI 
is used to create a superposition of two paths labelled 
|0) and 11). When the particle initially prepared in state 
|0) enters BSI, it transforms to |+) = (|0) + |I))/v^. 
Both paths are guided towards BS2 using mirrors. After 
the phase-shifter Z, the state of the particle is |+i) = 
(|0) + i|I))/\/2, and after BS2 it becomes {(I + i)|0) + 
(I — i)|I)}/2. Thus the particle has equal probability of 
reaching either of the detectors. The state |+) can also 
be written in terms of |±i), i.e., 

I \ 1 ~ E -v 1 ffi V -v /-I \ 

1+) = ^-i+*) + (1) 

We notice that the first component, namely |+i) trans¬ 
forms to |—) = (|0) — |I))/\/2 after Z, and then to |I) 
after BS2, and finally ends up in detector DI. Simi¬ 
larly, the second component, namely |—i) transforms to 
1+) = (|0) ffi |l))/v^ after Z, and then to |0) after BS2, 
and finally ends up in detector DO. In this sense, a mea¬ 
surement outcome of |0) (or |I)) amounts to postselecting 
\—i) (or l+i)) state just before the phase-shifter. 
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FIG. 1: (a) Three quantum particles entering a Mach-Zehnder 
interferometer consisting of two beam-splitters (BSl and BS2) 
and phase shifter (Z), and two particle-detectors DO and Dl. 
(b) Circuit for NMR investigation of QPHE. Hadamard gates 
perform the function of beam splitters, and Z-gate performs 
phase shift. Intermediate state information of the particle- 
qubits (Fi, F 2 , F 3 ) is encoded onto an ancilla qubit (H 4 ) 
using one of the controlled operations f/ 12 , f/ 13 , and t/ 23 . The 
ancilla qubit is measured at the end of the circuit. 


Let us now imagine the case of three particles in 
state 1000) entering the interferometer. After BSl, 
the state of the particles is described by the super¬ 
position \ 2 pa) = !+,+,+). The final state may col¬ 
lapse with equal probability to any one of the states 
{| 000 ), | 001 ), | 010 ), | 011 ), | 100 ), | 101 ), | 110 ), | 111 )}. 

The projectors 

P 12 = |0)(0| G) |0)(0| 1 + |1)(1| ^ 1 

Pi3 = |0)(0| G) 1 ^ |0)(0| + |1)(1| 1 ^ |1)(1| 

P 23 = 1 G) | 0 )( 0 | O | 0 )( 0 | + 1 G) | 1 ){ 1 | O | 1 ){ 1 |, ( 2 ) 

probe if any two of the particles are in the same state, 
i.e., 100) or |11). The expectation values of the pro¬ 
jectors give corresponding probabilities. Evaluating the 
expectation values for the state |+,+,+), we find that 
{P 12 ) = {P 23 ) = {P 13 ) = 1/2. Just after BSl, the prob¬ 
ability for any two particles being in the same path is 
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EIG. 2: The molecular structure of 3-bromo-2,4,5- 

trifluorobenzoic acid. The chemical shifts (diagonal elements) 
and effective coupling constants (off-diagonal elements) 
are shown in table. The NMR spectrum of ancilla is shown 

in Fig. I^a). 


therefore 1/2. 

We shall now consider only the cases wherein all the 
three particles reach the same detector, say DO (or Dl). 
Then the measurement outcome |000) (or |111)) is equiv¬ 
alent to postselecting the state |(/)i) = |—i, —i, —i) (or 
|0o) = |+L+L+^)) before the phase-shifter. 

The projection = Pjk\i^a) describes the com¬ 

ponent of lip a) corresponding to particles j and k being 
in the same path. Since, 

= (—L “'^1A 2 I+,+,+) 

(-i, -i, -i|Q, 0, +) + (-i, -i, -i|l, 1, +) 

2 

= 0, (3) 

we conclude that the postselected state |0o) has no com¬ 
ponent having particles 1 and 2 in the same path. Owing 
to the symmetry in the pre- {\ipa)) and post- (|0o)) se¬ 
lected states, the above conclusion can be extended to 
any pair of particles. This effect is interpreted as - if all 
the three particles have to reach the same detector DO, 
then no two particles can take the same path. Of course, 
similar interpretation can also be given for the case in 
which all the three particles reach the same detector Dl. 
This phenomenon which seems to violate the classical 
pigeonhole effect is called QPHE. 


3. NMR SIMULATION 

The 4-qubit quantum register consists of 3-bromo- 
2,4,5-trifluorobenzoic acid (see Fig. partially ori¬ 
ented in a liquid crystal N-(4-methoxybenzaldehyde)-4- 
butylanline (MBBA). The three ^^F spins and the 
spin are respectively used as the three particle qubits and 
the ancilla qubit. All the experiments were carried out 
on a 500 MHZ Bruker NMR spectrometer at an ambi¬ 
ent temperature of 298 K. The effective couplings in this 
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system are due to scalar interactions {Jij) as well as par¬ 
tially averaged dipolar interactions {Dij). Thus the in¬ 
ternal Hamiltonian for the system, under weak-coupling 
approximation, can be written as 

H = -27r ^ Vil^i + 27r y] Jljlzilzj, (4) 

i ij 

where Ui are the resonance frequencies, Jh = (Jij^2Dij) 
are the effective coupling constants, and I^i is the z- 
component of spin angular momentum operator of ith 
spin [21]. All the relevant Hamiltonian parameters like 
chemical shifts and effective coupling constants are tab¬ 
ulated in Fig. 

The first step in the experiment is the initialization of 
the 4-qubit system in a partial pseudopure state 

p( 0 ) = 1/16 -I- e| 000 )( 000 | (g) cr ^/2 ( 5 ) 
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FIG. 3: The ancilla NMR at various stages of QPHE sim¬ 
ulation, each obtained by a final 90° detection pulse. The 
simulated and experimental spectra are shown in the left and 
right columns respectively, (a) Thermal equilibrium state peq 
(the background baseline is due to the liquid crystal signal), 
(b) The partial pseudopure state p(0), and (c) after the entire 
MZI-circuit (without Uij) indicating the various combinations 
of detections. Spectra in (d-f) correspond to the complete 
QPHE circuit shown in Eig. mb) obtained with I/ 12 , U 13 , 
and U 23 respectively. The dashed boxes highlight the peaks 
corresponding to the postselected states as described in the 
text. 


EIG. 4: (a) Realizing Uij by a pair of GNOT gates, (b-d) 

NMR pulse sequences corresponding to G 1 NOT 4 , G 2 NOT 4 , 
and G 3 NOT 4 respectively. All the tt pulses (open rectangles) 
are about y axis and the phases of the 7r/2 pulses (filled rect¬ 
angles) are as indicated. The delays are set according to ti — 
1/(8J'4). 


where e ^ 10“^ is the purity factor. This state can be 
very easily realized by taking difference between two ini¬ 
tial states: (i) peq corresponding to thermal equilibrium 
state and (ii) pin obtained by inverting the populations 
of levels 10000) and |0001) using a low-power transition- 
selective Gaussian pulse of duration 80 ms [22] . The an¬ 
cilla spectrum corresponding to the initial state obtained 
by a 90° detection pulse (shown in Fig. [^(b)) clearly in¬ 
dicates the transition corresponding to the particle-state 
| 000 ). 

The next step is to apply the various elements of the 
Mach-Zehnder interferometer as shown in Fig. [^(b). All 
the unitary operations were realized by GRAPE optimal 
control technique [23]. The durations of the GRAPE 
pulses ranged from 400 ps to 700 ps, and the fidelities 
were better than 0.99 over a 10 % RF inhomogeneity 
range. 

The experimental spectrum after the complete MZI- 
circuit is shown in Fig. [^c). To probe the component 
of the wave-function prior to the phase gate, we insert 
operator Uij as shown in Fig. [^(b). The action of Uij is 
given by 

Uij = Pij (8) la + (1. - Pij) G) Xa, (6) 

where Pij are the projections as described in eqns. 1^, 
la are respectively the identity operators on system and 
ancilla, and Xa is the NOT operator on ancilla. Effec¬ 
tively, Uij conserves the ancilla if particles i and j are 
in the same path, but inverts it otherwise. Thus in the 
ancilla NMR spectrum, a positive ancilla peak indicates 
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Arrangement 

Possibility 

Classical 

Quantum 

DO 

Postselection 

Dl 

Postselection 

Any two in same 
container 

At least 
one pair 

None 

None 

Any two in different 
containers 

At most 
two pairs 

All three 
pairs 

All three 
pairs 

Any pair in 
container-0 
(container-1) 

Probable 

NO 

NO 

All three in same 

container 

Probable 

NO 

NO 

One of the particles 
in container -0 
(container-1) 

Probable 

YES (NO) 

NO (YES) 


TABLE I: The classical and quantum possibilities are tabu¬ 
lated for various arrangements of three particles in two con¬ 
tainers. The top two rows correspond to QPHE. 


two particles being in the same path, while a negative 
peak indicates them being in different paths. Each of 
the Uij operations is realized by a pair of CNOT gates 
as described in Fig. (a). The NMR pulse-sequences to 
generate each of the CNOT gates are shown in Fig. 
(b-d). 

The NMR spectra obtained after the complete MZI- 
circuit along with probing of particles ( 1 , 2 ) (by t/ 12 ), 
(1,3) (by t/ 13 ), and (2,3) (by I/ 23 ), are shown respec¬ 
tively on the right side of Figs. |^d-f). The correspond¬ 
ing simulated spectra are shown on the left side. The in¬ 
tensity variations in the experimental spectra are mainly 
due to the imperfections in executing CNOT operations 
resulting from RE field inhomogeneities, nonlinearities of 
the RE amplifiers, as well as decoherence. However, we 
observe an overall agreement of the experimental results 
with the quantum theoretical simulations. In particular. 


we focus on the NMR transition 000 (111) correspond¬ 
ing to postselection of all the three particles reaching DO 
(Dl). These transitions are highlighted by dashed lines 
in Fig. [^d-f)). All these transitions have negative in¬ 
tensities indicating that no two particles are in the same 
path. This is a clear demonstration of QPHE. 

For three particles and two containers, the various 
other arrangements and their classical as well as quan¬ 
tum possibilities are shown in Table 1. Here the first 
two arrangements, where there is a clear contradiction 
between classical and quantum, correspond to QPHE. 
For the other cases, while the possibilities are prob¬ 
abilistic in the classical regime, they are certain in 
the quantum regime. Of course, it is also possible to 
generalize QPHE to N particles in M < N containers [T]. 

4. CONCLUSIONS 

Quantum pigeonhole effect is yet another illustration of 
quantum systems displaying effects beyond the classical 
predictions. The nonclassical effect in this case emerges 
as a result of assignment of premeasurement states based 
on the knowledge of postmeasurement values. 

Here we provided the first experimental simulation of 
QPHE using an NMR quantum simulator. The quan¬ 
tum register consisted of three ^^F spins simulating three 
quantum particles whose intermediate state was probed 
by a ^H spin (ancilla). The experimental results agreed 
well with the quantum theoretical predictions. The suc¬ 
cessful demonstration of QPHE also illustrated good 
quantum control achieved on a four-qubit heteronuclear 
NMR system partially oriented in a liquid crystal. 
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